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PREFACE 

As  indicated  in  the  title,  this  paper  considers  the  problem  of 
the  Submarine  Commanding  Officer  who  must  determine  salvo  size 
and  spacing  in  making  a  torpedo  attack  on  a  target  vessel. 

The  paper  presupposes  a  knowledge  of  certain  concepts  of 
the  theory  of  probability.     The  terminology  and  notation  used  in 
discussing  xnathematical  probability  correspond  to  the  expressions 
in  general  usage  in  the  literature,    such  as,   "The  Theory  of 
Probability",  M.   E.   Munroe,    1951.     The  paper  is  offered  parti- 
cularly as  a  consideration  for  possible  classroom  discussion  by 
operations  analysis  students. 

The  i^riter  wishes  to  thank  Professor  Thomas  E.    Oberbeck 
of  the  U.    S.    Naval  Postgraduate  School  for  his  original  suggest- 
ions,  assistance,   encouragement  and  cooperation  in  the  prepara- 
tion of  this  paper. 
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GLOSSARY  OF  NOTATION 

A    s    track  angle 

B    a    target  true  bearing 

C    s    target  true  course 

D    s    domain  or  area  of  uncertainty  associated  with  position  of  MOT 

'*  S    ~    probability  density  functions 

h    z    hit  probability 

i    =    aim  point  designator;  i    =    1*2, ,n 

Ij    s    target  length 

a.  p.  8.     =    aim  point  sequence 

n    =    number  of  aim  points  in  an  aim  point  sequence 

m    =    number  of  hits  required  for  a  specified  degree  of  dannage 

s    nvimber  of  T/O's  fired  at  a  single  aim  point, 
p      s    minimum  acceptable  probability  of  a  hit 
R    ~    target  range 
S    s    target  speed 
T    =    time  on  station  in  days 
V  =    maximum  salvo  size 
N.    =    number  of  T/0*s  on  board  on  jth  day  of  patrol:  j  =  0, 1, ,T 

J 

N      =    number  of  T/O's  with  which  S/M  begins  patrol 
F.    =    expected  frequence  of  attack  on  jth  day  of  patrol 
B    -    angle  on  the  bow  (AOB) 
^  =    maximum  expenditure  of  T/O's  in  one  attack  on  jth  day  of  patrol 

J 


CHAPTER    1 
ELEMENTS  OF  THE  PROBLEM 
1,  1    Statement  of  the  Problem 

Each  time  a  submarine  commander  decides  to  make  an  attack 
he  must  make  a  complete  firing  phase  decision.     This  consists, 
among  other  things,  of  deciding: 

(1) 

(a)  How  xnany  hits  are  wanted  on  the  target. 

(b)  The  price  in  torpedoes  he  is  willing  to  pay  for  these  hits. 
Contained  also  in  the  firing  phase  decision  is: 

(2) 

(a)  Selection  of  the  final  inputs  into  the  problem. 

(b)  Selection  of  a  torpedo  spread  technique  to  insure  hitting 
the  target. 

(c)  Control  of  the  firing  of  the  torpedoes  (particularly  as 
to  time). 

This  paper  will  try  to  give  aid  to  (1)  by  an  analysis  of  (2) 
based  upon  statistical  inferences  concerning  2(a)  and  the  mathema- 
tics of  probability  concerning  2(bJ  and  2(c). 

The  basis  for  decision  will  be  derived  in  terms  of  probability 
of  hits  obtained  for  torpedoes  es^pend^d. 

The  handling  of  one  attack  is  then  to  be  correlated  with  the  con- 
duct of  a  complete  patrol  based  upon  what  is  desired  to  be  acabl^lished 
by  this  patrol,   as  specified  in  the  patrol  orders. 


1.2    Patrol  Objectives 

A  patrol  is  considered  as  taking  station  in  a  geographical  area 
under  one  of  two  alternatives  as  to  length  of  stay: 

1.  No  restrictions  as  to  time. 

2.  Stay  on  station  a  specified  number  of  days. 

Each  alternative  is  pursuant  to  patrol  orders.     The  first  envi- 
sages pursuing  every  opportunity  to  attack,   and  early  expenditure 
of  torpedoes  relieves  the  submarine  of  further  commitments  in  the 
area.     The  second  reflects  a  need  to  keep  the  submarine  in  the  area 
for  a  definite  additional  reason,   (e.g.    reconnaissance,    rescue,   or 
to  provide  a  constant  menace  to  the  enemy  over  a  specified  period 
of  time,  with  availability  of  relief,   say,   a  governing  factor). 

It  will  be  considered  that  submarine  patrol  orders  contain  the 
order  to  "Inflict  Maximum  Damage  on  the  Enemy"     which  may  be 
interpreted  to  specify  two  clearly  distinct  tasks: 

A.  "Insure"  kill  of  each  target  attacked, 
or 

B.  "Insure"  damiage  to  each  target  attacked. 

Here  the  word  "insure"  is  used  in  a  probabilistic  sense,   as  are  also 
the  words  "kill"  and  "dajnage". 

It  is  desired  that  the  submarine  commander  do  his  best  with  his 
opportunities  in  any  one  of  four  categories  in  which  his  orders  place 
him,  namely  Al,  A2,   Bl,   B2,   as  shown  in  the  following  table. 


1.     See  Theodore  Roscoe,   Submarine  Operations  in  World  War  II,    1949: 
Sub  Pac  Operational  Plan,   p.    233. 
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Kill  each  target  attacked.  Damage  each  target  attacked. 

Leave  when  torpedoes  Lieave  when  torpedoes 

expended.  expended. 


Space  the  kills.  Space  the  damages. 

Stay  on  station  a  specified  Stay  on  station  a  specified 

length  of  time.  length  of  time. 


In  a  type  2  patrol,   i.e.  ,  when  A2  or  B2  is  assigned,  note  that  it  is 
desired  that  the  commander  spreadj^.  his  effort  over  the  time  speci- 
fied so  that  his  ship  presents,   in  effect,   a  continuous  threat  to  the 
enemy.     [This  reasoning  should  also  leave  him  with  some  torpedo 
capacity  to  protect  himself  in  a  tight  spot  throughout  the  duration 
of  a  type  2  patrol.  ] 

It  will  be  seen  that  in  a  type  2  patrol  the  subznarine  Command- 
ing Officer  must  have  definite  knowledge  of  the  expected  frequency 
of  enemy  target  traffic  through  his  patrol  area  upon  which  to  base 
his  estimate  of  expected  frequency  of  attack.     This  implies  that  there 
will  be  sources  of  intelligence  from  which  he  can  estimate  this  target 
frequency  and  also  that  he  must  have  some  measure  of  his  own  search 
ability.     It  will  be  assumed  that  he  cank   from  previous  sources  of 
intelligence  of  whatever  origin,   (e.g.    intelligence  reports,   previous 
contacts,   previous  attacks,   etc.  )  make  such  an  estimate  of  the  ex- 
pected frequency  of  attack  on  any  given  day  of  his  patrol. 


1.  3    Firing  Phase  Objective 

For  uniformity,   we  define  the  firing  phase  as  the  period  from 
the  time  of  the  last  observation  (last  set  of  target  data  inputs)  until 
the  time  when  the  last  torpedo  of  the  salvo  crosses  the  target  track, 
real  or  extended. 

In  any  attack,  the  final  firing  set-up,   determined  by  whatever 
nneans,   is  put  into  the  torpedo  data  computer  at  a  given  instant. 
The  input  data  to  be  scrutinized  here  are:    True  bearing,   range, 
angle  on  the  bow(or  true  course)  and  speed  of  the  target.     True 
refers  to  the  compass  rose  only;  the  input  data  are  to  be  consider- 
ed as  a  set  of  random  variables  in  the  probability  sense.     The 
possible  errors  of  these  input  data  are  evaluated  to  predict  the 
size  and  probability  distribution  of  the  location  error  of  the  target, 
i.  e.   the  nature  of  an  area  (or  domain)  of  uncertainty,   say    D    , 
which  is  moving  and  expanding  with  elapsed  time,   and  into  which 
the  torpedo  spread  is  fired.     In  particular,   it  is  desired  to  know 
how  unit  mass,   representing  the  totality  of  probability  of  target 
position,   described  by  the  probability  density  function    f(D)  asso- 
ciated with    D    ,   is  distributed  over  this  area    D    at  a  certain 
elapsed  time    t    after  the  final  inputs  are  made. 

From  knowledge  of    D    and    f(D)    at  any  given  elapsed  timie    t  , 
it  is  desired  to  predict  the  aim  point  sequence,    spread  spacing,  and 
salvo  size  necessary  to  obtain;  (1)  at  least  one  hit  with  a  minimum 
^acceptable  probability    p      when  damage  is  wanted:    or,   (2)  at  least 
m    hits,  with  a  minimum  acceptable  probability    P       ,  when    m    hits. 


on  the  average,   are  needed  to  sink  the  target. 

The  Commanding  Officer  then  completes  the  firing  phase  deci- 
sion on  the  basis  of  the  information  cited  above.     If  a  type  Z  patrol 
is  being  conducted,   it  is  also  necessary  to  specify  the  salvo  size  for 
a  given  attack  that  is  compatible  with  the  requirements  of  the  length 
of  stay  on  station.     The  following  chapters  discuss  the  elements  of 
this  decision,   i.  e.    salvo  size  as  a  function  of  the  choice  of  one  of 
the  four  patrol  objectives  and  "how  good  the  set-up  is". 


CHAPTER    2 
ELEMENTS  OF  ONE- DIMENSIONAL  CONTINUOUS 
PROBABILITY  DISTRIBUTIONS 
2.  1    Scope 

Certain  fundajnentals  of  probability  will  be  discussed,  using 
simple  models,   in  order  to  develop  some  concepts,   definitions 
and  notation  for  later  use. 
2,  2    Target  Position 

Consider  a  target    T    of  length    L    on  course    C     .   For  conven- 
ience,    C       will  be  taken  as  the  positive  direction  on  the  axis  of  a 
o 

one -dimensional  coordinate  system.     Suppose  that,   at  the  instant    t, 
the  MOT  (middle  of  target)  is  known  to  be  on  the  x-axis,  but  its 
exact  position  is  uncertain.     The  position  of  the    MOT     ,   at  the  in- 
stant   t    ,   may  be  regarded  as  a  stochastic  variable  with  an  asso- 
ciated probability  density  function    f(x,;t).     The  function    f(x,  :t) 
will  be  termed  the  target  position  density  function  (at  time    t  ). 
In  the  examples  considered  below,   it  will  be  convenient  to  distin- 
guish between  the  cases  where    f(x.  ;t)    is  uniform  and  where 
f(x,  ;t)    is  normal. 

Hence,   when    f(x.;t)  is  uniform,    set    f(x.  ;t)    =    f  (x,  ;t)    where 

[eT-       ^°'"    ^o  -  7  ^^1    ^    ^  -^    ^o  +  7  ^^1 
(2.1)  yxj:t)     =      7       ^ 

I    0  elsewhere 

and    a      is  the  value  of    x.     denoting  the  mean  of  the  distribution. 
(See  Fig.    2a).     In  this  context,   target  location  error  is  taken  to  be    Ex. . 


When    f{xj;t)    is  normal,    set    fl^iJt)      =      f   (x,;t)    where 

i_         ,._._     /     1       /""r^o 


.    2 
X,  -a 


(2.2)  f(xi;t)     =     i—  Exp         -4  -4—^ 


where     -  oo  <    x,     <    <» 

^ain,     a     is  the  mean  of  the  distribution  but  in  this  context  the 

1 
target  location  error  is  taken  to  be    26'x,   . 

2.  3    Torpedo  Position 

Let    x,    represent  a  point  on  the  same  x-axis  where  a  hypothe- 
tical torpedo  (hereafter  called  T/O)  crosses  the  axis  at  the  instant 
t    ,     Then    h(x2;t)    ,  the  probability  of  a  hit  on  the  target  may  be 
computed  for  any  point    x^;  i.e. 


(2.3)  h(x2;t)      = 


The  function    h(x2;t)    is  shown  in  Fig.    2  b  for  the  case  where  the 
target  position  error    Ex.     exceeds     L.     In  this  connection,    set 
Ex,     =    M*L    where    M    is  a  constant,     M    >    1. 

Suppose  it  is  known  that  a  T /O  aimed  at  any  point  on  the  x-axis, 
say    x_    will  definitely  cross  the  x-axis  at  the  instant    t    ,  but  the 
exact  point  of  crossing  is  known  only  uncertainly  as  a  stochastic  vari- 
able with  the  associated  conditional  probability  density  function 
g(x2\;t|x2). 

Then,     P(x-;t)    ,  the  conditional  probability,  of  a  hit  on  the  target 
when  the  T  /O  is  aimied  at  the  p  oint    x-    is  given  by  the  formula 


1.     Some  attributes  of  the  normal  distribution  are  discussed  in 
Appendix  A 


(2.4)  P(x,;t)      =       I  h(x-,:t)  g(x,;t|x,)  dx,  . 


J' 

By  way  of  ilhistrating    P(x-;t)     ,  assume    g(x2;t|x,)    is  also  a 


uniform  probability  density  function  such  that: 


g(x2;t|x^)      =      \ 


li-  ^°^    ^3    -    I^^2   ^    ^2    ^    ^^3    ■*■    I  ^^2 

Ex^ 


0  elsewhere 

In  this  context,   the  quantity    Ex-,    will  be  termed  the  T /O  error. 
Assumie  also  that    f(x,;t)    is  given  by    (2.  1).     Finally,   let 
Exj    >    Ex^   ^    L  ,   i.  e.,     Ex.     =    M»L    ,     Ex^    =    NtL    ,  with 
M    >    N    >    1    .     The  graph  of    P(x-;t)    under  the  assumptions  of 
this  paragraph  is  shown  in  Fig.    2  c. 

Whenever  the  instant    t    is  not  essential  to  the  discussion,   it 
will  be  convenient  to  write 

f(xj;t)       =       f(x^) 

Hx^rt)     =       Mx^) 

P(x3;t)      =       P(x3) 
etc. 

Since  the  primary  interest  of  this  chapter  is  the  probability 
density  fnnction  associated  with  target  position,   it  will  be  assumed 
that  the  T/O  error,     Ex^     ,   is  zero.     Under  this  assumption, 

(2.  5)        j      x^  =        x^ 

.3)         =        hCx^) 


2.4     Computation  of  hit  probabilities. 


An  hypothetical  example  is  inserted  at  this  point,  primarily 
as  a  vehicle  to  explain  the  notation  used  in  the  papero 
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It  is  convenient  to  formulate  firftt  of  all  a  system  for  designating 
the  crossing  points  of  T/O's  on  the  k-axis.     By  assumption  [see  (2.  5)  J 
the  crossing  point  of  the  T/O  is  the  aim  point  of  the  T/O.     The  choice 
of  aim  points «  and  the  probable  coniequences  thereof  is  a  basic 
objective  of  this  paper.     A  definite  choice  of  eq\ially  spaced  aim 
points  will  be  called  an  aim  point  sequence  (a.  p.  s.  ).     Inasmuch  as 
it  will  be  assumed  throughout  this  pa. per  that  all  T/O's  cross  the 
axis  at  the  same  instant    t  ,  the  a.  p.  s.    is  not  concerned  with  the 
order  of  firing;  it  specifies  only  the  number  of  aim  points,  where 
aimed,  and  the  spacing  between  aim  points.     Thus,  let  an  a.  p.  s. 
be  defined  as  follows: 

a      s  point  marking  th^  mean  of   f(x,) 

d      ^1/2  distance  between  adjacent  aim  points  along  the  axis 

n      sl: number  of  aim  points  in  the  a.  p.  s. 

See  Fig.    2  d. 
Then,  according  as    n    is  Odd  or  even  the  aim  points  in  an  a.  p.  d.  are: 
n    odd  n    even 

^1    =  ^o  ^1    =  Vi^ 

a,      sr.a      +    2d  a,     =    a^    +    d 

2  o   —  2  o 

a.,      =    a      T    2d  a,      =    a_    +    3d 

3  o  .   3  o  — 

a.      =    a     T    3d 
4  o 

etc. 

Attention  is  turned  fsrst  to  the  case  where  a  single  T/O  is 
placed  at  each  aim  point  of  an  a.  p.  s.  ,  with  spacing    2d    =    L  , 
a  so-called  simple  salvo.     Later,  the  case  of  placing  more  th|m  one 
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T/O  at  each  aim  point  of  an  a.  p.  s.   will  be  discussed. 

For  the  case  of  a  simple  salvo,    let    P.    denote  the  probability 

of  a  hit  associated  with  the  i  th  aim  point  of  an  a.  p.  s.  ,     i  =  1 ,  .  .  .    ,  n. 

n 

Let                       P,(n)    =        U  P.       . 

^                   i=l  ^ 

i.  e.  ,     P,  <n)    is  the  total  probability  of  one  hit  for  a  sinmple  salvo 

of    n    aim  points. 

By  way  of  illustration.   Table  IJgives  values  of    P.     and    P,(n) 

for  simple  salvos  in  two  cases.     The  first  case,   designated  in  the 

table  by    f       ,   is  for    f{x, )    uniform  with    Ex^    =    4L    [cf.    formula 

(2.  1)]  and  the  second  case,  designated  by    f       ,   is  for    f(x. ) 

normal  with   (Tx,     =    L    [cf.   formula  (2.2)] 

Table    I 


n  = 

1 

n  = 

2 

n  = 

3 

n  = 

4* 

a.  p. 

f 
u 

f 
n 

f 
u 

f 
n 

f 
u 

f 
n 

f 
u 

f 
n 

^1 

.250 

.383 

.250 

.341 

.250 

.383 

.  250 

.341 

^2 

.250 

.341 

.  250 

.242 

.  250 

.341 

^3 

.250 

.242 

.  250 

.  136 

^4 

.250 

.  136 

Pi 

.  250 

.383 

.500 

.  682 

.  750 

.867 

1.000 

.954 

In  the  case  of  simple  salvos,  with    2d    =    L,  note  that 
there  are  no  gaps  nor  overlaps  in  the  spread  spacing.     This 
insures  far  any    n  ,  a  maximum  value  of    P,     in  the  case  of  sym- 
metrical target  position  density  functions,   such  as    f  (x,)    and 

f„(Xj). 
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Two  possible  criteria  for  the  choice  of    n    are  as  follows: 

(a)    Let    p      3    minimum  acceptable  probability  of  a  hit  with 

one  T/O.     In  this  case,     n    =    maximum    i    for  which    p.^  p     , 

'^i       '^o 

(b)    Let    P      minimum  acceptable  total  probability  of  one  hit 
for  a.  p.  s.    selected.     In  this  case,     n    =    minimum    i    for  which 

By  way  of  illustration,    suppose    p      is  fixed  at  0.  2;  then  from 

Table  1,  for    f       say,  it  can  be  seen  that,     n    =    4  would  not  be 
n         ' 

used,  but  instead    n    =    3  would  be  used.     If,   however,   a  high 
P      was  desired,   say    0.  9  ,  then    n    =    4    would  necessarily  be 
used.     The  choice  between  (a)  and  (b)  is  a  matter  of  preference. 
Hereafter,   the  choice  of    n    will  be  given  in  terms  of  both 
criteria. 

The  possibility  of  more  than  one  hit  with  a  salvo  to  T/O's 
arises  if 

(1)  more  than  one  T/O  is  fire^  at  each  aim  point  of  an 
a,  p.  s.   with  the  distance  between  adjacent  aim  points 
exactly  equal  to  target  length,   i.  e.  ,     2d    =    L  ; 
or  (2)   one  T/O  is  placed  at  each  aimi  point  of  an  a.  p.  s. 

with    the  distance  between  adjacent  aim  points  less 
than  target  length,   i.  e.  ,     2d    <    L  ; 
or  (3)    more  than  one  T/O  is  fired  at  each  aim  point  and 

2d    is  less  than    L  . 
These  possibilities  will  be  considered  in  the  order  just  enumer- 
ated and  will  be  referred  to  as  Cases  1,   2,  and  3,   respectively. 
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These  cases  will  be  considered  under  the  simplifying  asstunptions 
that 

(I)  each  T/O  crosses  the  x-axis  exactly  at  its  intended  aim  point. 

(II)  all  T/O's  of  a  salvo  may  be  considered  to  cross  the  x-axis 
at  virtiially  the  s^me  instant    t    ,   even  if  more  than  one 

T  /O  is  fired  at  any  aim  point. 

In  considering  Cases  1,   2,  and  3,  let    ^(n,m,2d)  denote  the 
probability  of  obtaining  two  hits  with  a  salvo  of  size    nxm    in  which 
m    T/O's  are  placed  at  each  of  the    n    aim  points  of  an  a.  p.  s. 
with  spacing    2d    between  aim  points. 

Under  assumptions  (1)  and  (U)  above,   it  follows  that  in  Case 

1  there  are    m    hits  if  there  is  one  hit  and  hence  that 

P_(n,m,L)    =    P,(n)  ;  that  is»  the  value  of    P     (n,m,  L.) 
m  1  m 

is  identical  with  the  value  of    P,(n)    for  a  simple  salvo.     For 

example,   referring  to  Table  1,   if  3  T/O's  are  aimed  at    a.     , 

then  three  hits  would  be  obtained  with  probabilities    .  250  and 

.  383  ,   respectively,  for  the  target  position  density  functions 

f      and    f      of  Table  1.     Similarly,  for  a  salvo  of  6  T/O's,   it  is 
u  n  ' 

found  that    P2(3, 2,  L)  has  the  values    .750    and    .867,   respec- 
tively, for  the  target  position  density  functions    f      and    f      of 
Table  1. 

Turning  to  Case  2,  the  probability  of  more  than  one  hit  depend? 
in  part  upon  the  amount  of  "overlap"^  between  successive  aim  points 
of  the  a.  p.  s.     To  illustrate,  consider  a  salvo  of  two  T/O's  fired 
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in  an  a.  p.  s.   of  two  aim  points  with  spacing    2d    =    "T  L.     Then  for  a 
target  position  density  function    f  (x,)    with   (T        =    L    (as  in  Table  I) 

Pj(2,  1,^L)    =    .547 

P^lZ.l.-jL)    »    .197 
Finally,  in  Case  3,   it  is  clear  that  if    m    T/O's  are  fired  at 
each  aim  point  of  the  a.  p.  s.   with    2d    <    L  ,  there  is  a  possibility 
of    m  ,   2m,   etc.   hits.     In  this  case  it  is  possible  to  calculate 
P.     (n,m,  2d)    where    j    is  a  positive  integer  having  values  1,   2,  etc. 

•I 

up  to  some  maximum  value  depending  upon    2d.     For  a  given  salvo 
of    nxm    T/O's,   it  is  important  to  note  that  whenever    n  ;>  1    and 
2d    <    L    ,  a  penalty  is  paid,   namely,   reduction  in  the  value  of 
P     (n,mi,2d)    to  achieve  the  probabilities    P.     (n,m,2d)    where 

3^  1. 

The  spacing    2d    =    L    has  been  measured  along  the  x-axis, 

coincident  with  the  assumed  target  course    C     .     Consider  the 

T/O's  on  any  track    T      other  than  normal  to  the  x-axis.     Let 

o 

A   denote  the  track  angle  between    T      and    C      (measured  as  in 

»  o  o    ^ 

Figa   2  e),  where  for  simplicity  adjacent  tracks  are  considered 
parallel).     To  intercept  an  interval    L    on  the  x-axis,  the  tracks 
must  be  spaced    JL  sin  A    apart.     L  sin  A    is  called  the  effective 
target  length.     For  the  one -dimensional  case,  it  is  clear  that  the 
probability  of  a  hit  associated  with  a  given  aim  point  does  not 
depend  upon  the  track  angle.     To  avoid  consideration  of  target 
width,  very  sharp  track  angles  will  not  be  considered  here 
(roughly,     A  <  15      or    A  <  165       is  not  considered). 
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f  (x,:t) 
u     k 


1/E: 


o       Z         1 


^0+7^^1 


Figure  2a:    f^(xj;t) 


Mx^;!) 


o 
(M-l)L 

ML 
(M+1)L 


Figure  2  b:    h(x^;t)    where     Ex,     =    ML  ,   M>1 


P(x3;t) 


Figure  2  c:     P(x3;t)    where     Ex       =    ML,     Ex^    =    NL  ,   M>N>1 
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^n<-l) 


n  =  1  - 

n=  2  . 
n=  3. 
n  =  4. 


■^ 


Aim  Point 
Sequence 
(a.  p.  8. ) 


Figure  2  d:    The  a.p.  s.    depends  upon    n    ,   the  number  of  aim  points, 
and  the  probability  density  associated  with  target  position. 


Figure  2  e:    Aina  point  spacing  along  the  normal  to  T/O  track. 
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CHAPTER    3 
ELEMENTARY  PROBABILITY  DENSITY  FUNCTIONS  FOR  A 
STOCHASTIC  VARIABLE  WHICH  DEPENDS  LINEARLY  ON  TIME. 

Suppose  that  at  a  certain  instant  of  titne    t    ,  measurements 
are  made  of  the  position  and  speed  of  a  point  which  is  moving 
along  a  straight  line  at  constant  speed.     Let  the  position  and  speed 
xneasurements  both  be  subject  to  error  and  suppose  that  these  mea^ 
surements  may  be  regarded  as  observations  of  two  independent 
stochastic  variables.     Then  it  is  frequently  necessary  to  "predict" 
the  position  of  the  point  at  some  later  time    t  +  T.     In  this  con- 
text, "predict"  means  to  determine  the  probability  density  func- 
tion associated  with  the  position  of  the  point  at  time    t  +  T. 
The  future  position  of  the  point  is  a  stochastic  variable  which 
depends  linearly  on  time.     If    X    denotes  the  future  position  of  the 
point  at  time    t  +  T  and  if    x    denotes  the  position  of  the  point  at" 
time    t    ,  then 

X    =    X    +    sT 
where    s    is  the  speed. 

This  chapter  discusses  the  probability  density  function  asso- 
ciated with    X    when  the  probability  density  functions  of    x    and 
8    are  (1)  both  uniform  and  (2)  both  normal. 

3,1    Probability  Density  Functions  of    x    and    s    Both  Uniform. 
Given;  (a)    A  uniformi  probability  density   f ,    ■    f-(x;t)    in 

the  interval    (x      -    7  £x  ,  x      +   -j  ^^^    associated  with  the 
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xneasurexnent  of  the  position  of  the  point  on  the  line  at  time    t    . 

(b)    A  uniform  probability  density   f,    "    fJi^lt)  in  the 
interval    (s^    -    -y  E     ,8      -f   •«- £   )    associated  with  the  speed  of 

O  ^0  O  M        0 

the  point  on  the  line.     When  the  position  of  the  point  at  time    t    is 
known  exactly  and  when  time  is  considered  as  being  measured 
without  error,  then  the  position  of  the  point  after  an  elapsed 
time  y  is  uncertain  and  the  associated  probability  density 
function,  here  denoted  by    f^    ^    f2(8T;t)    ,  is  uniform  in  the 
interval    (s^T  -    "jEsT,     s^T     +   ■jEsT') 

Clearly  at    'y  s    0  ,    f,(x;t)    completely  describes  the  pro- 
bability density  associated  with  the  point*  with  mean  at    x 
and  density    f.     s    -= —  in  the  interval    (x      -    7  Ex,  x      +    7  Ex) 
and    0    elsewhere.     This  would  be  the  case  at  any  instant    t 
when  a  measurement  of    x    is  made.     If  no  further  measurement 
of    X    is  available,  but  the  position  of  the  point  at  time 
t.     =    t  +  ^  is  required,  then  it  is  necessary  to  consider  the  new- 
stochastic  variable 

X    =    X  +  sr 

The  probability  density  function    f^    of  the  new  stochastic  variable 

1 
is  obtained  by  the  convolution  of  f;.    and    f^  . 

Thus,  fj    =    f3(X;tj)    =    fj(X;t)    *    f2(X;t) 

»     f^(X-  8T)f2<«T)rd8 

=   ^f^CX  -  x)  f ^(x)  d^ 
For  the  simple  case  here,  the  distribution  diagram  for    f^    is 
represented  in  Fig.    3  a. 
1.     Cf:    M.   E.   Munroe,   Theory  of  Probability,    1951,  Section  38. 
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It  is  geometrically  evident  that: 

(1)  The  mean  of    X    is  moving,  and  at  time    t.    is  at    x      -f    u/y  • 

i  o  o ' 

(2)  The  values  of    X    are  restricted  to  lie  within  the  interval 

t^o  +  ^"^  -  4^^  +  4^-'^'  ^o  +  •o'^  +  4ex  +  ^E»r)J 

which  lengthens  withf  , 

(3)  The  area  of  the  trapezoid  representing  the  probability  density 
function   f^    over  this  interval  miust  always  equal  one  (the  total-^ 
ity  of  probability).     Consequently,  the  ordinate  values,  which 
represent  probability  density,   must  necessarily  be  reduced  as 

T'    increases. 

3.  2    Probability  Density  Functions  of   x    and    s    Both  Normal. 

Consider  now  the  case  wbiere  each  of  the  density  functions    f^  , 
{y  ,  f .    is  not  uniform,  but  normal.     In  accordance  with  the  nota- 
tion in  Appendix    A 

^1      '     ^^^'^o'   k'^^ 

Then  for    any  time    t.    =    t  +  T »  the  stochastic  variable 
X    =    X  +  s  T  has  a  normal  probability  density  function  which, 
by  the  addition  theorem  for  normal  distributions  miay  be  repre- 
sented in  the  notation  of  Appendix    A 
as  fj      5     f3(X:X^,  (T^it^) 

Where       X       =      x      +    s    T 
o  o  o  ' 

(^2     ,    ^  2    ^   ^2^2 

X  x  s 

The  consequent  movement  of  the  normal  variate    X    ,  the 
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expansion  of    (T      and  the  thinning  out  of   f  ^    is  also  clear  in  this 


X 

case  as  shown  in  Fig.   3  b. 


1^ 


f^w 


Es 


1/Es 


i 


f3(x) 


^  2a_J 

2h-2a 

LZa  _^ 

/ 

l/2b 
if 

\ 

X  +8   r 
o     o 


2b  =  Max[Ex.,E8T 
2a  =  Min[Ex      Est 


Fig.    3a:    f,    =    f     *f      when    f.     and    f^    are 

both  uniform 


/,(Xj) 


■sr 


fjCx) 


(  rxj)^  +  (<r8r)^ 


Fig.   3  b;    i. 


X  +   8     T 
O        O 


f,     *    £2    when    f.     and    f^    are 


both  normal 
20 


CHAPTER   4 

THE  AREA  OF  UNCERTAINTY  OF  TA|IGET  POSITION  AND  ITS 

ROLE  IN  THE  SUBMARINE  FIRING  PHASE  DECISION. 

4,  I    Preliminary  Considerations  and  Notation 

Let  the  independent  stochastic  variables    (R«B)    describe  the 
position  of  a  point  in  a  plane;  e.g.  ,  an  observation  of  the  range    R 
and  the  true  bearing    B    of  a  target.     Let  the  uncertainty  of  the 
exact  value  of    R    (at  time    t    )    be  associated  with  the  uniform 
probability  density  function 


fi(R;t) 


for    R^    -    J  ER   S    R   :i    R^    +   -j  ER 


0  elsewhere 

Let    B    be  associated  with  the  uniform  probability  density 
function 


fjCBU)    . 


^  for    B^    -   ^EB:S    B  i   B^    +   |  EB 


I 


0  elsewhere 

For  definite ness,  consider    B    measured  in  radians,   say  from  a 
stationary  submarine  whose  position  will  be  considered  a  fixed 
datum  point  of  the  problem  in  any  coordinate  system  desired. 
The  Joint  probability  density  function  associated  with  a  pair  of 
values    (R»B)    at    t    is  then 

f3(R.B:t)    =    fj(R;t)  f2(B;t)  . 
Thus,    f^    is  positive  over  the  annular  sector  shown  in  Fig.   4  a 
as  the  domain    D,    =    D,(t)  . 
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Suppose  now  that  no  ctlatter  what  point  in   D.    is  selected,  t&ere 
is  a  vector    V-,  uncertain  as  to  direction  and  speed,  that  can  send 
this  point    (R,B:t)    into  a  new  position,  say   (P«/^t,T)    after  an 
elapsed  time  'y    o     As  heretofore,  consider  ^    as  stop-watch  time 
measured  from  the  instant    t    without  error.     As  shown  in  Fig.  4  b 
consider  the  angle    B     that  'O*  makes  with    B    to  be  an  indepen- 
dent stochastic  variable,  with  associated  probability  density 

function 

^    ^  '"    \    "   1"'^  ^    "  ^    "o    '^  7 

0  else>flrhere 

Consider  the  speed    S    ,  i.  e. ,  the  magnitude  of    '\r  ,  to  be  an  in- 
dependent stochastic  variable,  with  an  associated  probability 
density  function 


_    1  ««.«.«.     1 

f5(S;t) 

0  elsewhere 


for    S      -    iES    <    S  ^    S      +    iES 
o         2  "^  o  2 


Then  at  the  end  of  a  specified  time    'T^  ,  the  probability  density 
associated  with  the  magnitude  of  the  vector    l/'T'  due  to  the  variate 

"^   '"  %ST.t)    =    (^         '"    ^°    -    T^ST  i    Sr  <    S^    +  I  EST 

1    0  elsewhere 

Evidently,   for  each  possible  origin    (R,  B;t)    in    D^  ,  the  vector  of 

uncertainty  "^T^   has  a  domain  of  termination,   say 

D,     »    D,  (R,B;"V',  Y  I*)    consisting  of  the  annular  sector    D, 

shown  in  Fig.   4  b.     And  with    D,     is  an  associated  probability  density 

function 

f^(s.e:r:t)  =  f^ceu)  fgCST.t) 


22 


3ut  "Vf  may  originate  irom  any  pojnt  in    D-     and  further  "XrT  may 
assume  any  one  of  the  spt  of  values  adlowed  by    D,     .      It  may  be  ct>nsl- 
de red  that    D^    "maps"    D^    into  a  new  domain    D_  «  D-(P,^  ;t,T) 
where   v     and  p    are  new  stochastic  variables.     Functionally,     D_ 
may  be  represented  as 

(4.1)       D^    =:    D^  [R.f^(R);B,f2(B):e.f4(e);S.f5(S):r:t] 
The  region    D-    may  be  obtained  as  the  result  of  a  simple  geome- 
trical construction  and  is  an  eight-sided  figure  (see  Fig.   4  c). 
The  boundary  of  the  area    D-    is  comprised  of  the  four  sides  of 
D-    and  the  four  sides  of    D.   .     With    D-    is  associated  a  positive 
valued  probability  density  function    f_    =    iji(^ »P  itty)    .       The 
mathematical  problem  is  to  deterniine  exactly  how  unit  mass, 
representing  the  totality  of  probability,   is  spread  over    D_;  f_ 
denotes  the  description  of  this  mass  distribution  over    D_    . 

No  attempt  to  develop  the  analytic  theory  for  the  evalua- 
tion of  the  joint  density  function  will  be  made  here.     Clearly,   the 
probability  density    f^    of  a  set  of  values  of  the  new  variates 
{P,fi)    at  time    (t+T)    is  zero  at  the  boundary  edges  of    D_    and 
rises  continuously,  with  possibly  a  different  slope  from  each 
edge  to  a  plateau  of  maximum  probability  density,   say 

f_    =    max  f- 
It  is  evident  that  there  exists  no  simple  way  of  describing  the 
infinite  variety  of  configurations  of  D_    and  the  associated  proba- 
bility density  function    f „    for    D_    and    f_    depend  upon  the  rela- 
tive magnitudes  of  all  of  the  parameters  indicated  in  formula  (4.  1). 


23 


For  example,   Fig.   4  d  shows  the  regipn    D      which  results  when    D^  , 
the  area  of  uncertainty  of  target  positipn  at  time    t    ,   is  propogated 
with  no  speed  error  and  a  discrete  course  error. 

As  a  point  pi  reference  in    D„  ,   define    A       ^o  be  the  center 
of  gravity  of  the  mass  distribution  represented  by    f_  ,     The  coordln-  • 
ate s  of    A      will  be  designated    (P    ,^     )    or    (P    ,0    ;t,T)    when  the 
time  parameters  are  essential.     When    \    and    ^     are  independent 
stochastic  variables  then   \        =  r      and    ^       -  p     where    \     and  ^ 
denote  mean  values  of    P    and   p     respectively. 

Alternatively,   if    f_    has  a  maximum  value  at  a  unique  point 

A      of    D-    ,   then  this  point  may  be  selected  as  a  point  of  refer-    - 

ence  of    D^    and  in  this  case  the  coordinates  of    A      will  be 
»  o 

designated    {P    t^   )    where 

VP*.^*)    =    Maxf^Cp,/?)  . 

(e,/S)€D^ 

In  the  familiar  special  case  where    f-    is  known  (or  approximated) 
in  the  form 


(4.^)       f;(x,y)    =   3^^^^,  .  Exp. 


(x-x^)         (y-Yo' 


«   y  2<r2  26-^ 

then    A      =    (x   ,y   ).     (See  Appendix  A,   Fig.   A-2).     In  this  special 
case  the  iso-density  contours  defined  by    f _    =  £       where    £     is  a 
constant  are  concentric  ellipses. 

Proceeding  frcm  the  formulation  of  D_  and  f_  with  the 
sinnplest  of  assumptions  concerning  the  component  stochastic 
variates    R,B,&,S,     it  would  be  desirable  to  turn  to  the  case 
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where  the  component  errors  are  independently  and  normally  distri- 
buted.    Tlje  reader  will  recognize  the  inherent  xnatKcmatical  difff- 
culties  which  are  beyond  the  scope  of  thi*  paper. 
4.  2    The  Submarine  Fire  Control  Problem 

The  concepts  and  notation  of  Section  4.  1  will  now  be  related  to 
the  submiarine  fire  control  problem.     Thus,  let  '*(R,B;t)    denote 
the  actual  range  and  bearing  of  MOT  at  the  time  of  final  observa-' 
tion  of  the  target  (instant  of  last  inputs  into  Position  Keeper  of 
TDC  before  firing).     R      and    B       are  the  mean  or  expected  >|alues 
of    R    and  .B    at  time    t    with  respect  to  the  distributions    f.(R;t) 
and    f2(B;t)    ,   respectively.     The  functions    f,     and    f,    are 
assumed  to  be  known.     It  will  be  assumed  further  that  the  values 

of    R      and  B       are  identical  with  the  final  observations  of    R 
o  o 

and  B      (inputs  into  TDC).     Similarly,     Q    and    S    denote  the  AOB 

(angle  on  the  bow)  and  the  speed  of  the  target  at  time    t  .     6     and 

S      are  the  mean  or  expected  values  of    6^    and    S    at  time    t    with 
o 

respect  to  the  functions    iA^lt)    and    fc(S;t)    respectively.     Again, 

the  functions    i.     and    fc    are  assumed  to  be  known  and  the  values 
4  5 

of    B      and    S      are  assumed  to  be  identical  with  the  final  obser- 
o  o 

vations  of    6"    and    S    (inputs  into  TDC).     Then  the  region    D- 
may  be  called  the  area  of  uncertainty  of  MOT  at  the  tiraie  of 
last  observation  of  the  target.     The  region    D_    then  represents 
where  the  MOT  may  be  found  at  time    t+T   (y units  of  time  after 
the  last  observation)  when  there  are  uncertainties  in  target 
range,  bearing,   course,   and  speed.     The  region    D_    will  be 
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termed  the  area,  of  uncertainty  of  MOT  at    t+  T'  . 

The  intervi^l   'Y  may  be  considered  as  the  suxn  of  two  intervals 

ri  and  Xz  '       '^  =  T-i  +  r^ 

The  interval    t',,   represents  the  "dead-time"  between  the  last 
observation  of  the  target  and  the  instant  of  firing;  the  interval 
Tl     ,  which  is  a  function  of  H^,     )s  the  tiine  required  for  a  "perfect" 
T/O  to  reach    a^    S    (?Q,^^:t,T)    when  the  T/O  is  fired  at  the 
instant    t+fi  •     Briefly,    7^2    ^^  T/O  running  time  and 

The  submarine  is  assumed  to  have  a  TDC  which  can  predict 
1^    within  the  interval    y,     ,  i.  e,  ,   T/O  running  time  is  calculated 
before  the  T/O  must  be  fired. 

To  this  point,   the  direction  of  the  ta^rget's  znotion  has  been 
discussed  in  terms  of    6    ,  the  AOB .     It  will  be  convenient  to 
speak  now  of  the  target's  course    C     ,  where,   of  course, 
C    =    C(9,B)    and  in  particular    C      =    C(6    ,B   ).     Inasmuch  as 
C      is  the  target's  course  as  it  appears  in  the  TDC  in  accordance 
with  previous  assumptions,  it  is  convenient  to  consider  the  aim 
points  of  the  eum  point  sequence  (see  Ch,    2)  of  the  T/O  salvo  to 
be  along  the  course  line    C       extending  through  the  point    a 

Unlike  the  simple  one -dimensional  case  discussed  in  Chap- 
ter 2,   there  is  in  the  two-dimensional  problem  no  simple  way 
to  calculate  the  probability  of  a  hit    P.    associated  with  the  ith 
aim  point  of  the  a.  p.  s,   when    f„    is  generated  by    f.,  f,,  f^, 
and    fc    .     To  discuss  the  two-dimensional  case  further,   it  is 
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convenient  to  use  the  simplest  of  continuous  probability  density 
functions  as  a  model  for    f-  . 


To  simplify  them.^consider 

for 

^o  -  i  ^y  -  y  '  ^o  "^  i  ^y 


ICTV         ^^'^    ^o-    f=^    ^    *^    ^o+    f=^ 


f^    =    f(x.y)    = 

'  0  elsewhejre 

With  this  simplification,     D-{P  ,S  ;t,7^)    is  a  rectangle  which  will  be 

designated    DCt+f)  .     See  Fig.   4  e. 

Recall  that    D    represents  the  area  of  uncertainty  of  position 

of  MOT  at  time    t+  T'.     The  time  interval  is  computed  for  the 

collision  instant  of  a  "perfect"  T/O  with  the  point    a     ,     Let    C      » 

•  o  o 

the  "most  likely"  target  course-   through    a    ,be  identical  with  the 

X-axis  of  coordinates.     Let  the  collision  instant    t+'T'  =    t      .     To 

o 

simplify  further,   assume  that  all  T/O's  have  infinite  speed  through 
D(t   )    ,   and  that  all  T/O's  of  a  salvo  pass  through    D(t   )    at  the 
same  instant    t       .     The  width  of  the  target  will  be  ignored.     For 
T/O's  traversing  paths  parallel  to  the  y-axis,  the  probability  for 
a  hit  at  a  given  aim  point    X^    along    C      would  be 


h(x2;t^)    = 


^2    +    jL 


^  I     f(x,  y)  dydx 


; 


^2    "    I  ^  ^-*» 


which  reduces  to  the  result  (2.  2)  for  the  one -dimensional  case  of 
Chapter  2,     For  track  angles  other  than  90     this  is  not  the  case. 
Referring  to  Fig.   4  f ,   it  will  be  seen  that  tfee  probability  of  a  hit 
for  such  a  track  is  simply  the  volume  under    f(x,y)    "swept  out" 
by  the  T/O  "path"  of  width    L  sin  A    . 
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Consider  now  the  assumption  that  the  T/O  has  infinite  speed  fn 
the  light  of  the  simple  example  just  discussed.     Taking  into  account 
the  finite  speed  of  the  T/O,   a  very  precise  calculation  of  probabi- 
lity of  a  hit  would  require  determination  of  the  earliest  possible 
instant,   say    (^q-^i)    after  firing  at  which  the  T/O  enters  a 
region    D(t   -  S, )    in  which    f(t   -  S. )    has  a  positive  value  and  the 
latest  possible  instant,   say    (t  +  S'yf   at  which  the  T/O  leaves  a 
region    D(t  +5^)    in  which    fU+S",)    has  a  positive  value.     When 
f(t   -  6". )    and    Ut'i'Sy)    differ  sufficiently  over  the  "path  swept 
out"  by  the  T/O,  there  may  be  an  appreciable  error  in  assum- 
ing that  the  T/O  passes  with  infinite  speed  through  a  region 
D(t   )    and  sweeps  out  a  path  under  the  probability  density  func- 
tion  f(t   )  .     Moreover,  the  statemejits  above  that  the  T/O 
enters    D(t^-S"j)    and  leaves    D(t^+5^)    imply  that  these  regions 

have  certain  boundaries  which  will  not  be  true  when    f    has 
positive  values  over  the  entire  plane  as,  for  example,  when  at 
least  one  of  the  component  probability  density  functions    f. ,  f,  * 
f.  ,    or    fe    is  normal.     Under  these  circumstances  it  maybe 
convenient  to  consider  a  region,   say    D(t   -5",)    bounded  by  an 
C-valued  iso-density  contour  of    f(t  -o.)  (cf.   remarks  follow- 
ing formula  (3.  2))  which  the  T/O  enters  at    t   -o,     and  another 
region,   say    D(t  +  6,)     ,  bounded  by  an    ^-valued  iso-density 

contour  of    f(t  +5".}    which  the  T/O  leaves  at    t  +5",    and  deter- 
o      c,  o      c 

mine  whether    fft   -0.)    and    f(t  +«,)    differ  "significantly" 

OX  G  £t 

over  the  path  swept  out  by  the  T/O,     When  these  two  functions 
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do  not  differ  significantly,  then  the  probability  of  a  hit  associatecF 
with  an  aim  point  ntiay  be  visxialized  as  the  ii|olutne  under    f(t  :^) 
swept  out  by  the  path  of  the  T/O  through    D(t   :£)  . 

In  actual  practice,  the  probability  density  functions    f, ,  1^$ 
i.  t     and    fc    associated  with  target  range,  bearing,  course,  and 
speed  inputs  may  not  be  known.     However,  it  seems  evident  from 
the  analysis  of  this  paper  that  these  functions  are  essential  opera- 
tional data  which  are  necessary  for  the  submarine  fire  control 
problem  if  there  is  to  be  a  rational  basis  for  T/O  salvos.     It  is 
likely  that  all  four  functions  are  normal  density  functions  with 
standard  deviations  which  maybe  estimated  by  operational  trials. 

Appendix  B  gives  a  brief  discussion  of  the  numerical  approxi' 
mation  of    f-    by  high  speed  digital  computers  when    f .   ,  f ,  ,  f . 

and    fc    are  given. 
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D2(R.B;t)  =  a  b 

a  b  =  c  d  =  ER 

a  d    s  (R     -  7  £R)EB 


b  c  «  (R     +  J  ER)EB 


f3(R.B;t)    « 


Figure  4  a:    The  annular  sector    D,(t)    and  the  associated 
probability  density  function    f  ,(R,  B;t) 
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D,     =    a  b  c  d_ 

P    =    (a    -    yE8)r    Ee 

Q    =    (8     +    ^E8)r    EG 


Figure  4  b:    The  annular  sector    D,  (R,  B;v',T;t) 
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Figure  4  c:    The  region    D^ 
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D_    s    abed 


Figure  4  e:     The  case  where    D_    is  a  rectangle 

f(D)    '<      ^^3cE 
<L        0 


E3EEV       ^^   °7 


elsewhere 


Figure  4  f:    The  effect  of  T /O  track  angle  on  the  area 
of    D_    which  is  "swept  out". 
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CHAPTER    5 
GENERAL  l^EMARKS  ON  TO|lPEDO  SALVOS 

5.  1    Torpedo  Lethality 

The  number  of  T/0'»  required*  on  the  average,  to  destroy    a 
target,  depends  on  two  primary  factors: 

(1)  The  accuracy  of  the  T/O,  i.  e.  the  probability  the  T/O 
will  hit  the  target,  and 

(2)  The  destructiveneis  of  the  T/O,  i.e.  the  probability 
that  the  target  is  destroyed  if  the  T/O  hits  it. 

Here  (1)  is  reduced  to  the  problem  of  finding  the    P.'s    , 
thence    n    . 

For  (2)  in  many  ca«es  a  single  hit  ii|  not  enough  to  destroy  a 
target.     Even  for  merchant  ships  a  singly  T/O  hit  is  not  usually 
enough  to  sink  it,  and  heavy  combatant  s^ips  are  designed  to 
withstand  many  T/O  hits.     To  tjceat  such  cases  exactly,  one  should 
determine  the  probabilities,   say,     d,   ,   d2  ,  d^  ,  etc.   of  sinking 
the  ship  if  the  ship  is  hit  1,   2,   3,   etc.  .  times.     Then  if,  for  a 
given  method  of  aiming  T/O's,  the  probabilities  of  1,   2,   3,   etc. 
hits  are    P,   ,  P,  «   P^  »   etc,  the  probability  of  sinking  the  ship 
would  be    d,P,  +  d2P2  "*"  ^\^%  "*" »    where  the  d's  would  repre- 
sent damage  coefficients. 

For  this  paper  it  is  simply  assumed  that  there  is  a  number    m 

for  each  type  ship  such  that,  on  the  average,     m    hits  of  that  type  will 

n 
fifink  thft: Siai)  with  probability    P      =    P(ip,n,  d)=   Tl  P. 

°  i=l  '        -^ 
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with    d    =    ly  L  sin  A      ',;  the  spacing  along    the  normal  to  T/O  tracks 
for  the  a.  p.  s.  described  in  Chapter  2. 

5.  2    Time  Intervals  in  Salvo  Fire 

r  '■' II. 

Recall  that,  in  firing  t;orpedoes  fronri  i  submarine,  the  firing 
is  not  a  true  spread  in  the  sense  that  all  t/O's  are  launched  at 
once  to  provide  a  simultaneous  menace  to  tk^  enemy.     T/O's  are  ' 
fired  singly  and  consecutively  with  a  prescribed  time  delaiy  to 
avoid  countermining. 

Let    A*V  s    time  delay  (an  interval  of  ieveral  sees. ).  Consi- 
der a  spread  of    nxm    T/O's,  let    -^  •    1  .     L«t  7^     be  initial 
make  ready  time  after  last  inputs  to  TDC  at  time    t  .     The  mean 
T/O  would  go  at  time    'y   ♦   -ii^i— AT   after  time    t    ,  on  an 
expected  run  time    ^T?    8*y»     Let  mean  ^.    =  T*      +    -—r —  ^^  '» 
then  mean  T^  =Tj    +^2      '^^^   *o    ~    *    **"  '^* 

The  mean  track    To  ,  commenced  at  time    t  +  7i    may  be 
considered  as  the  datum  line  about  which  the  spread,  or  aim 
point  sequence  (a.  p.  s. )  is  placed. 
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CHAPTER    6 
THE  FIRING  PHASE  DECISION 

6.  1    The  Cardinal  Rule 

The  number  of  aim  points    n    in  an  a.  p.  s.   miay  be  regarded 
as  depending  upon  the  function    f(t   ;^)    introduced  in  Section  4.  2. 
Hence    n    is  a  function  of     'y^ , 

Therefore,  as  a  oghrdinal  rule  to  minimize  the  spreading  and 
thinning  out  of  probability  density,   avoid  protracted  firings: 

FIRE  AS  SOON  AS  POSSIBLE  AFTER  DECISION  TO 
COMMENCE  FIRING  PHASE. 
6.  2    The  General  Procedure 

Most  of  the  notation  used  here  has  sJready  been  discussed, 
otherwise  it  is  defined  in  the  glossa;ry.     Four  sets  of  procedures 
are  outlined,  following  the  dictates  of  the  four  objectives  pro- 
posed in  Chapter  1.     In  each  case    n  ='maximum    i    for  which 

P.   ^    p        and    d    =    -y  L  sin  A    along  the  normal  to  T/O  tracks 

n 
in  the  aim  point  sequence.     Then    P      =    ^     P. 

i=l 
I.     Fire  Salvos  to  metximize  probability  of  at  least  one  hit 

(i.e.   case  Bl).     Do  not  fire  T/O's  at  aim  points  for  which  asso- 
ciated probability  of  hit  is  less  than    p     .     As  nearly  as  possible, 
all  T/O's  are  to  be  expended  in  not  more  than    T    days;  obtain 
extra  hits  as  indicated  by  frequency  of  attack.  f 

(a)    If  ]2  ^    "     ^i^®  o'^e  T/O  at  each  of  first  ^j  aim  points 
(a.  p. )  (cf.   Section  2.  3) 

'    -i     :  ■        •  -:        I        ri--^  ■'.' 

\  4  *"  "  .   .  .  w  .         •  -.*  --      .    p- 
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(b)  If    n   <  J]    lookatC'      then  if : 

} 

(i)  C   <Cj     fi'^e     *A  T/O's  at  each  of  first    n    a. p.    and 
one  T/O  at  each  of  first  ^   -    /xn     a.  p. 

(ii)    n  ^Ci    "^G'    ^^"^^     u  T/O' s  at  each  of  first    n    a.  p. 
(iii)  ^  ^'    '^    T^   <  ^    fi*"©  one  T/O  at  each  of  first    n    a.  p. 

(c)  When  2j  '      "    ^^   request  permission  to  reduce  value  of   p 

(d)  Return  to  base  when  all  T/O's  are  expended*  but  do  not 
remain  on  station  more  than    T    days. 

II.     Fire  T/O's  to  maximize  probability  of    mi    hits  on  targets 

requiring    m    hits  to  be  sunk  with  probability    Po    (i.e.  «  case  Al)". 

Do  not  fire  at  aim  points  for  which  associated  probability  of 

hit  is  less  than    p      . 
o 

As  nearly  as  possible,  all  T/O's  are  to  be  expended  in  not 
more  than  T  days;  obtain  extra  hits  as  indicated  by  frequency 
of  attack. 

(a)  If    mn  ^2    »  ^"^    "^    T/O's  at  each  of  first   \)    a.  p. 
(Optional  to  fireC-    ml)   T/O's  at    (i)*  l)th    a.  p.). 

(b)  If    mn    <  YL  »  determine  ^3  '       then  if: 

(i)   £7  ^  Sj     ^^^^  yu.  T/O's  at  each  of  first    n    a. p. 
(ii)    mn<  C  •     <  C    ^^^^    >*•'    T/O's  at  each  of  first 

n    a.  p. 
(iii)  2'.     <    mn   <2      fire    m    T/O's  at  each  of 
first    n    a.  p. 

(c)  When  ^'.       a    2^    request  permission  to  increase 

•I 

value  of    Po    and /or  to  reduce  value  of    p     . 

o 
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(d)    Return  to  base  when  all  T/0*«  are  expended^  but  do  not 
remain  on  station  mor^  than    T.     days, 

III. Remain  on  station  for  approximately  T.  days,  condxicting 
as  many  attacks  as  possible.  Ffre  salvos  to  maximige  probability 
of  at  least  one  hit  (i.  e. «  case  B^ 

Do  not  fire  T/O's  at  ainm  points  for  which  associated  probabi- 
lity is  less  than  p  .  As  nearly  as  possible,  expend  all  T/O's  by 
end  of  patrol;  obtain  extra  hits  as  indicated  by  frequency  of  attack. 

(a)  If    S,    ^   n   fire  one  T/O  at  each  of  first  £?.    a.  p. 

(b)  If    n   <  £J^     fire  yuL^  T/O's  at  each  of  first    n    a. p. 
and  one  T/O  at  each  of  first   2      -  Mn    a.  p.    . 

(^)    When  y7  \     *    ^  IZ    'c<iu«*t  P«"^*il<>'^  to  '^cluce  value 

of   F|i- 

(d)    Return  to  base  when  all  T/O's  are  expended,  but  make 

this  as  close  as  possible  to    T.    days. 

JY.     Remain  on  station  for  approximately    T.    days,  conduct- 
ing as  many  attacks  as  possible.     Fjre  salvos  to  maximize  probabi- 
lity of    m    hits  on  targets  requiring    m    hits  to  sink  with  probability 
P   (i.  e.,    case  A2). 

Do  not  fire  T/O's  at  aim  points  for  which  associated  probabi- 
lity of  hits  is  less  than   p     .     As  nearly  as  possible  expend  all 
T/O's  by  the  end  of  the  patrol;  obtain  extra  hits  as  indicated  by 
frequency  of  attack. 

(a)    If  2Ji    -^    nan    fire    m    T/O's  at  each  of  first  l),    a.  p. 
(Optional  to  fire  TJl  "    "*^1    T/O's  at    (l)  ^    +  1)^   a.p. ). 
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(b)  If   mn   <   2j.    fir«yiLT/0's  at  each  of  first    n    a.  p.  (Optional 
to  fire  J3j    -  /*i^   T/0'»  at    (n  •!•  1)  th_  a.p. ) 

(c)  When  ^3  *    '    ^]^  request  permission  to  increase  value  of 

P      and /or  reduce  value  of    p_  . 
o  o 

(d)  Return  to  base  when  all  T/Q's  are  expended,  but  xnake  this 
as  close  as  possible  to    T.    days. 
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Appendix    A 
The  Normal  Distribution 

A-1    The  One-Qiinensional  case. 

The  normal  distribution  is  one  of  the  most  significant  mathema- 
tical models,  from  both  the  theoretical  and  practical  standpoints, 
upon  which  the  techniques  of  applied  statistics  are  founded. 

Using  as  an  example,   say,   a  set  of    n    measurements,   with 

errors    x.(i  =1,2, ^)  ,     some  underlying  assumptions  of  the 

distribution  are; 

(a)  In  performing  the  set  of  measurements,   not  adl  the  errors 
X.    are  equally  likely  to  occur, 

(b)  In  general,   large  errors  are  less  likely  to  occur  than 
small  ones. 

(c)  If  there  is  a  best  value    m    ,  then  positive  and  negative 
errors  about  this  value  are  equally  likely  to  occur, 

(d)  If  (c)  is  true,   then    m    is  the  arithmietic  average  of  the 
individual  measurements. 

The  density  of  this  distribution  is  given  by 

Thus,   a  normal  distribution  is  completely  determined  by  specify- 
ing its  mean    m    ,   and  standard  deviation    u     i,  e,     f(x;m,(r). 
The  function  is  shown  in  Fig,    A- 1  for  several  values  of    <5~   , 
Changing    m    merely  shifts  the  curves  to  the  right  or  left  without 
changing  their  shape. 
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By  the  linear  transformation    y    =  ~-4 —    ,     f(x)    is  put  in    stan- 
dard form,  with    mi    =    o    and     0^  =    1  .     The  density  function  is  then 

f(x:o,l)    = Exp  (   -    yx     ) 

It  is  in  this  form  that  one  usually  finds  the  integral  tabulated  for 

any  value  of    x    . 

A- 2    The  Two -Dimensional  Case. 

A  simiple  case  of  two  independent,   normally  distributed  varia- 
bles will  be  discussed  here.     Given    f,(x;x   ,  (T   )    and    fo(y;y   tC^'S 

1         o       X  L  '   'o       y 

such  that    f(x,y)    =    f,(x)  i-Ay)    .     Let  the  point    (x   ,  y  )    =    a 
The  joint  density  is  given  by: 


^  X     y  ' 
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(i.e.   unit  mass  is  spread  over  a  symmetrical  mound).   Fig.   A-2, 

Curves  formed  by  a  vertical  slice  parallel  to  an  axis  are  of 

the  same  shape  as  the  normal  distribution  curve. 

The  curve  forined  by  a  horizontal  slice  _is  an  ellipse 

2  /        X       2 


f(x.y)    =; 


.x-a^   >^        ^      /y-a. 


^x 


=     c     (a  constant) 


The  equation  represents  an  ellipse  with  center  at    a      semi- 
axes  of  length    (c(S~  ,  c  (T  )    and  axis  parallel  to    (x,  y)    axels.      An 
ellipse  of  this  nature  is  called  an  iso-density  contour  ellipse. 

The  probability  that    (x,  y)    lies  inside  the  ellipse  equalis  P    when 

2       *A/  2  2 

c       =  /X^        ,     f    =    2    (where    c       is  found  using  the    x,  y    values  and 

P    is  the  probability  in  per  cent  for  the  corresponding  fractile  of 

the  Chi- squared  distribution  with  2  degrees  of  freedom). 
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To  define  a  bound  for  the  area    D    covered  by    f(x,y),   any 

contour  ellipse    C    maybe  selected,  e.g.   noting  that: 

2 


y 


\/^ 


90 


=     4.61       (f=2) 


then    c       =4.61  c    =    2.  15 

let    F     =    90%    contour  ellipse  with  semi-axes    (2.  ISCT  ,2.  15^  ) 

'^  X  y 

i.e.   for  £"    =    90%    the  ellipse  is    4.  3(r     by    4.3^  ,  centered  at    a   , 

•^  X       '  V  o 


Figure  A-1:    Normal  Frequency  Curves. 
m=0;  er=  0.4»  1.0,2.  5 
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Figure    A  -  2:    Bivariate  Normal  Distribution 


A-4 


Appendix    B 
Proposal  for  Preparation  of  Histograms  Used  in 
Salvo  Firing 

A  proposal  for  using  an  electronic  digital  computer  for  deter- 
mining   D    and    F(D)    suitable  for  use  as  reference  diagrams. 

First,  prepare    statistically,   the  histograms  representing 
R,     f(R);  B,   f(B):  e,  f(6)    where    6    =    0(0,  B);    and    S,  f(S). 
Thus  for    R,  we  might  have  the  histogram  shown  in  the  following 
figure. 


^1 


A  A3        A4       A5         A^       A^ 


^1         ^2        ^3        ^4       ^5        ^6       ^7 


For  each  component  the  rnachine  is  instructed  to  pick  a  random 
number  (by  a  s^lb- routine  method)  say    0    ^    r    i    1    and  this  num- 
ber instructs  the  machine  from  which  step  of  the  cumulative  graph 
to  pick  the  value  of  R.     e.  g.  ,   0    Si    r  Jk   A,     picks    R,,     A,    <    r   ^ 
A.     +    Aj    picks    Rp,   etc.     In  this  manner,   a  value  of  each  compon- 
ent is  picked  and  the  machine  computes  the  simple  trigonometry 
involved  to  get  the    X,    Y    coordinates,    say  from  the  datum  point 
for  a  given     ^, 
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These  are  compared  to  the  grid  values,    sorted,   and  added  to  the 
proper  cell  as  a  simple  count  of. the  occurrence 
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^3 

^4 
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^(^^4^2) 

^3 

^4 

^5 
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boxes 


D  and  f(D)  for  various  set  of  components  and  'T^'s  could 
be  make  up  to  guide  the  approach  officer  in  the  selection  of  the 
p.'s    and    n    . 

Not  only  would  one  like  histograms  based  on  various  mean 
values  to  use  in  the  computation,  but  also  on  the  method  of 
measurement,    such  as  range  histograms  at  various  ranges 
taken  by:    radar,  optical,    sonar,   etc:    Bearing:    optical,  hydro- 
phone,  sonar;  Speed;  screw  count;  tracking,   etc. 

If  each  histogram  was  made  up  in  standard  steps  the  compu- 
ting routine  could  be  used  over  and  over  again  to  construct    D's 
and    f(D)    for  various  concibinations  by  just  changing  the  values 
in  the  cumulative  histogram  cells. 

Note  that  the  computer  would  produce  discrete  mounds  that 
would  be  as  good  as  the  histogram  data  from  which  they  were 
prepared. 
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No  mention  of  what  the  component  error  data  actually  looks 
like  is  made  here,  because  of  the  non-classified  nature  of  the 
paper. 
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